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SEPARATING INEQUALITIES FOR NONNEGATIVE POLYNOMIALS 
THAT ARE NOT SUMS OF SQUARES 



O " SADIK ILIMAN AND TIMO DE WOLFF 

<N' 

Abstract. Ternary sextics and quaternary quartics are the smallest cases where there 
exist nonnegative polynomials that are not sums of squares (SOS). A complete classifi- 
cation of the difference between these cones was given by G. Blekherman via analyzing 
the extreme rays of the corresponding dual cones. However, an exact computational 
approach in order to build separating extreme rays for nonnegative polynomials that 
are not sums of squares is a widely open problem. We provide a method substantially 
simplifying this computation for certain classes of polynomials on the boundary of the 
PSD cones. In particular, our method yields separating extreme rays for every nonneg- 



■ ative ternary sextic with at least seven zeros. As an application to further instances, we 



compute a rational certificate proving that the Motzkin polynomial is not SOS. 



(N ■ l. Introduction 

>' 

We consider real polynomials in the vector space of all homogeneous polynomials in n 
variables of degree d, denoted by H n d . For every p e H n d we denote its real projective 
variety as V(p). Let P n ^ C H n d be the cone of all nonnegative polynomials in n variables 
of degree d. 

Inside H n ^di there are two full dimensional convex cones of special interest, the cone of 
nonnegative polynomials and the cone of sums of squares (for a general background about 
nonnegative polynomials and sums of squares see e.g. [TOj [111 HSl 120]; for some metric 
and convexity properties of these cones see [2]). 



£ i Pn,2d := {/ e H nM : f{x) > 0, for all x G W 1 } 



■>n,2d 



{/ e P n ,2d ■ f = fi for some f E . 



The investigation of the relationship between the cone of nonnegative polynomials and 
the cone of sums of squares began in the seminal work of Hilbert when he showed that 
the cone of nonnegative polynomials coincides with the cone of sums of squares exactly 
in the cases of bivariate forms (n = 2), quadratic forms (2d = 2) and ternary quartics 
(n = 3, 2d = 4) ([T3]). 

The Motzkin polynomial m(x, y, z) = x 4 y 2 + x 2 y 4 — 3x 2 y 2 z 2 + z 6 was the first explicitly 
known example for a nonnegative polynomial which is not a sum of squares. Most proofs 
for this fact are based on term by term inspections (see e.g. [171 20J ) . In near past other 
proofs were found, e.g. using representation theory (see [8]). 
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In |3J Blekherman showed that for fixed dimension 2d > 4 there are significantly more 
nonnegative polynomials than sums of squares as n tends to infinity. However, the question 
of precisely when nonnegative polynomials begin to significantly overtake sums of squares 
is much less understood. In the smallest cases where there exist nonnegative polynomials 
which are not sums of squares ((n, 2d) = (3, 6), (4, 4)) the general conjecture is that these 
two cones do not differ very much. This conjecture is supported by the following two 
facts: Firstly, the maximal dimensional difference between exposed faces of the cone of 
nonnegative polynomials and sums of squares is one (see [4]). Secondly, all extreme rays 
of the dual sums of squares cone Eg 6 have rank one or rank seven (see [5]). 

Recently, in [6] it is shown that except the discriminant there is a unique component 
of the algebraic boundary of S 36 with degree 83200 which indicates the complicated 
structure of the SOS cone. But still, the geometry and the relationship between these two 
cones in the smallest cases are less understood. 

In the smallest cases (n, 2d) = (3, 6) and (n, 2d) = (4, 4) Blekherman showed that 
it is precisely the Cayley-Bacharach relation that prevents sums of squares from filling 
out the cone of nonnegative polynomials. More precisely, in |5j it is shown that every 
separating extreme ray in the dual SOS cone for a given nonnegative polynomial that 
is not a sum of squares depends on a 9-point configuration for (n, 2d) = (3, 6) resp. a 
8-point configuration for (n, 2d) = (4, 4) coming from the intersection of two qubic resp. 
three quadric polynomials. Furthermore, given an appropriate 9-point (resp. 8-point) 
configuration, one can write down an extreme ray of the dual SOS cone (see Theorems 
12.61 and 12.7ft corresponding to faces of maximal dimension of the SOS cone. 

A central problem in this area is how to obtain the separating inequalities efficiently 
This can always be done in a numerical way (see Section |2T2~|) . but is widely open for exact 
methods currently. Hence, finding constructive methods for computing these inequalities 
is one main research issue. Blekherman's result does not provide an efficient symbolic way 
to obtain a proper 9-point (resp. 8-point) configuration to solve this problem (see 12.21 for 
further details). 

The key idea of this article is to construct a proper 9-point (resp. 8-point) configuration 
out of a given initial set of points. Specifically, we investigate nonnegative polynomials p 
which lie on the boundary of the cones P^q and (which cover most of the explicitly 
known nonnegative polynomials that are not SOS, see [3]). Our main result, Theorem 
provides a sufficient condition for using k zeros of p as a subset of a 9-point (resp. 
8-point) configuration. The idea is to fill up the set of k zeros with 9 — k (resp. 8 — k) 
points such that a genericity and a quadratic condition based on the Cayley-Bacharach 
relation holds (note that k < 10 for p G P-s^ \ S 3j6 and p G \ £ 4j4 ; see [6j [9]). 
Given these conditions, which are computationally easily checkable (at least for (n, 2d) = 
(3, 6)), we can construct a separating extreme ray immediately. This method reduces the 
complexity of constructing separating extreme rays via symbolic computation significantly. 
Furthermore, it yields rational certificates for rational point configurations and even for 
rational varieties V(p) C Q 3 in case of (n, 2d) = (3, 6). 
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We show that for p G ^3,6 \ £3,6 and k = 7 (resp. p G \ £4,4 and k = 6) almost 
every 9-point (resp. 8-point) configuration containing the seven (resp. six) zeros leads to 
a certificate for a nonnegative polynomial p to be not SOS. 

In Section [2] we review some curve theoretical issues as, e.g., the Cayley-Bacharach 
relation and present Blekherman's results on the SOS cones S 3 6 and £4,4. In Section [3] 
we state and prove our main Theorems 13.11 and 13.31 for ternary sextics and quaternary 
quartics and discuss exactness and rationality of our methods. Section @] deals with the 
special case of polynomials with exactly seven zeros. We show that in this case our 
method generically yields a separating extreme ray (Theorem 14. ip . Finally we discuss the 
difficulties of dropping zeros in our method in Section [5J via applying it to the Motzkin 
polynomial and constructing a rational extreme ray certificate for it. To the best of our 
knowledge, our certificate is the first separating rational extreme ray certificate for the 
Motzkin polynomial. Furthermore, based on the Motzkin polynomial we investigate some 
geometric aspects of the set of appropriate point configurations for our method (see Figure 
[TJ). In the appendix we discuss an example of the seven point case and an example of 
a point configuration for the Motzkin polynomial where our method does not yield a 
certificate. 



2. Preliminaries 

2.1. Curve theoretical background. We recall some classical results from algebraic 
geometry. From now on we consider every investigated polynomial to be homogeneous. 
We start with the Cayley-Bacharach relation. It exists in various formulations (see |12j); 
we use the ones given in [3]. 

Lemma 2.1. Let pi,p2 G C[xi, £2, £3] be two ternary cubics intersecting transversely in 
nine projective points 71, ... , 79. Let t>i, . . . , v 9 be affine representatives ofji. Then there 
is a unique linear relation on the values of any ternary cubic on Vi 

9 

(2.1) /] u jf( v j) = f or Q tt ternary cubics f 

3=1 

with nonzero Uj G C. Furthermore, if (12. ip is satisfied, then the following genericity 
condition holds 

(2.2) no four of the Vi lie on a line and no seven on a quadric. 

Analogously, we have for quaternary quadrics 

Lemma 2.2. Let p\,p2,P3 G C[x\, . . . ,£4] be three quaternary quadrics intersecting trans- 
versely in eight projective points 71, . . . , 73. Let vi, . . . , v$ be affine representatives of 7$. 
Then there is a unique linear relation on the values of any quaternary quadric on Vi 

8 

(2.3) u jfi v j) = f or a M quaternary quadrics f 
3=1 
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with nonzero Uj G C. Furthermore, if ( 12. 3 p is satisfied, then the following genericity 
condition holds 

(2.4) no five of the Vi lie on a plane (i.e. projective linear 2-space). 

For the genericity condition ( 12. 4 p see e.g. [THj. Note that, if all points Vj are real, then 
all Cayley-Bacharach coefficients Uj are real, too (see e.g. [5J Lemma 4.1]) and can be 
computed by solving a system of linear equations with the coefficients of forms in i7 3 3 
resp. H^2 as variables. 

Each of the conditions (I2.2p and (12. 4p can be checked easily by investigating the minors 
of the matrix given by the vectors Vj. 

For the description of separating extreme rays of Eg 6 yielded by Blekherman's Theorem 
!2.6l one needs to investigate 9-point configurations given by intersecting two ternary cubics. 
These ternary cubics need to be coprime. The following lemma shows that this generically 
is the case (see [2Tj). 

Lemma 2.3. Suppose A := {v\, . . . ,vs} is a set of eight distinct points in M. 3 , no four 
on a line and no seven on a quadric and let fi, f'2 be a basis for the vector space of all 
homogeneous cubics with projective variety affinely represented by A. Then fi and $2 are 
relatively prime. 

This lemma yields that one can apply Bezout's theorem in order to compute a ninth 
intersection point Vg of fi and f 2 . However, v Q might not be different from Vi, . . . ,v 8 (i.e. 
the intersection multiplicity might be greater than 1). But, again, generically this will 
not be the the following lemma shows (see [18]). 

Lemma 2.4. Let fi, f'2 be two homogeneous polynomials in n variables (with n > 2) 
of degree d and generic coefficients. The discriminant A(fi, f 2 ) vanishes if and only if 
fi{x) = f2{x) = has a singular solution. The set of polynomials for which this is the 
case is a hypersurface. 

In Section 0] we will investigate the special case of polynomials p G \ £3,6 with 
exactly seven zeros. In this context we use the following lemma (see [2l]). 

Lemma 2.5. Suppose A is a set of seven distinct points in M. 3 , no four on a line and no 
seven on a quadric with basis fi, f2, f3 for the vector space of homogeneous cubics with 
projective variety affinely represented by A. Then fi, f 2 , / 3 have no common zeros outside 
of A. 

2.2. Blekherman's results. In [5J Blekherman was able to fully characterize the extreme 
rays of the dual SOS cones E3 6 and £4 4 , basically via using the Cayley-Bacharach relation. 
We recall his main result, concentrating on ternary sextics. 

Theorem 2.6. Suppose that a ternary sextic p is nonnegative but not a sum of squares. 
Then there exist two real cubics q\,q2 G -£^3 intersecting transversely in 9 projective points 
71, . . . , 79 which yield a certificate forp G P 3i6 \£ 3i6 . More precisely, let v\, . . . , v g be affine 
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representatives 0/71, ... ,79. Then there exists a linear functional I : H^^ — > K given by 

9 

i(f) = J>/K) 



/or some a ; 6 C such that l(p) < and /(£ 36 ) > 0. Furthermore, at most two of the 
points 7i are complex. 

This theorem works the same way for p G P^ \ £4,4 (see [5]). 

Recall that for every I G S3 6 there is a corresponding quadratic for Q/ defined by 
Qi : P 33 — >• R, / !->■ /(/ 2 ) (see e.g. [SJ US])- One defines the rank of a linear functional 
/ G Sg 6 by rank(Z) := rank(Q^). In [5 J it is shown that every extreme ray of £ 36 , which 
does not correspond to point evaluation (i.e., a rank 1 form), is given by a rank 7 form 
which comes from a 9-point evaluation. 

Theorem 2.7. Suppose I spans an extreme ray o/Sgg which does not correspond to 
a point evaluation. Let Wi be the kernel of the corresponding quadratic form Qi and 
suppose gi, 02 G Wi intersect transversely in 9 real projective points 71, . . . , 79 with affine 
representatives v\, . . . ,Vg such that the unique Cayley-Bacharach relation is given by 

U 1 f(v 1 ) + '--+Ugf(Vg) = for f G Hs,3- 

Then Qi can be uniquely written as 

Qi(f) = a 1 f(v 1 ) 2 + --- + a 9 f(v 9 ) 2 

with exactly one single negative coefficient a& and the rest of the strictly positive and 
ak is given by 



(2.5) a k 



21? _| i_ Ha _ ?i 

ai ag ak 



Furthermore, any such form is extreme in £3 6 . 

Again, an analogue version for p G \ E 4 4 holds (see [5]). 



Suppose p G Ps^ \ S 3 6 and we want to construct a separating extreme ray I for p 
using the upper theorem. Therefore, we need to find two coprime ternary cubics 01,02 
intersecting in 9 points. But, q±, 02 need to be contained in the kernel Wi of the quadratic 
form Qi corresponding to /. Hence, one already needs to know / to determine q±, q2- 

This problem can be avoided by choosing a 9-point configuration A = {v±, . . . ,Vg} 
coming from an intersection of some real ternary cubics qi,q2- So, a separating extreme 
ray I is obtained by finding an appropriate a = (a\, . . . , ag) satisfying (12. 5p with respect 
to A such that l a (p) < 0. The question whether such an a exists is unclear a priori though 
it can be answered by quantifier elimination methods (see e.g. [U [7]). But, to the best of 
our knowledge, no methods are known to compute an appropriate a in a symbolic, exact 
way efficiently 
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However, one can solve this problem numerically. Let p G P3^ \ S 3 6 be a ternary sextic 
and r G int (£ 36 ) (e.g. r = x 6 + y & + z 6 or r = (x 2 + y 2 + z 2 ) 3 ). Consider the following 
semidefinite optimization problem. 

min A such that p + Ar G So K 

AeR d ' D 

For A minimal the polynomial p + Ar is strictly positive and lies on the boundary of 
S 3j6 . Hence, p + Ar is a sum of exactly three squares s 2 + s| + s| (see 0). 

The polynomials s\, s 2 , S3 have no common zeros and an appropriate linear combination 
of two of these polynomials can be used as q% and q 2 in Blekherman's theorem. Of course, 
the computation of the corresponding nine intersection points will be difficult and not 
exact, too. Furthermore, getting "nice" values (e.g. a rational minimal A) depends also 
highly on the choice of the polynomial r G int (E 3 6 ) . It is not clear how to choose r in 
dependence of p. 

In the case p G \ £4,4 this approach works the same way. For A minimal the 
polynomial p + Ar is a sum of exactly four squares p + Ar = s 2 + s 2 + s§ + s\. Three of 
these four Sj have a common zero (see [5]). 

3. A Certificate for boundary polynomials not to be SOS 

Our approach to construct a separating extreme ray for a given boundary polynomial 
p G dPsfi \ £3,6 is to investigate certain point sets A := {vi, . . . , v 9} containing the variety 
V(p), satisfying the genericity condition (12. 2p and for which we can certify that there are 
coprime polynomials qi, qi G H 3 3 with V(qi) H V{q?) = A. 

Note that, if we talk about zeros of homogeneous polynomials in this and the following 
sections, then we always consider their affine representatives with slight abuse of notation. 

The easiest case is when p has at least eight zeros V\,...,v% (satisfying (12. 2p ). Lemma 
12.31 provides the existence of coprime qi, q 2 vanishing on Vi, . . . , v 8 and thus a ninth point 
Vg is given by Bezout's Theorem. For a generic set of zeros Vi, . . . ,vg the corresponding 
coprime polynomials qx,q% have generic coefficients and hence, due to Lemma \2A\ we 
have v 9 {v±, . . . , v s} generically. Thus, A := {v 1, . . . , v 9} satisfies f 12 . 2 j) generically. This 
yields a certificate I immediately since for any choice of ai, . . . , a§ we obtain an ag < by 
(12. 5 p and hence 

9 

Kp) = ^2 a jP( v j) = a 9P(v 9 ) < 0. 

J'=l 

In the following we generalize this idea to any number of zeros between one and eight. 
We choose the zeros v±, . . . , v & of a polynomial p G dP^^ \ £3,6 as a subset of the nine 
intersection points A := {v\, . . . , vq} of two coprime ternary cubics. We provide a symbolic 
method based on genericity conditions which yields a separating extreme ray if one finds 
a (9 — fc)-point configuration satisfying some quadratic relation. Specifically, the following 
theorem holds. 

Theorem 3.1. Let p G dP%§ such that p is not SOS. Let A := {vi, . . . ,Vg} C M 3 such 
that the genericity condition (12.21) holds and V(p) = {vi, . . . ,Vk} with 1 < k < 8 (i.e.: the 
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projectivization of A is the intersection of two coprime polynomials qi,q2 G ^3,3/ Then 



there exists a certificate l a : H 36 — > R, / t-t Ylj=i a jf( v j)' 

) G M tyit/i 

respect to A for p to be not SOS , if the following inequality holds 

(3.1) (u 2 k+1 + ■ ■ ■ + u 2 8 )(p(v k+1 ) + ■ ■ ■ + p(v 8 )) < ulp(v 9 ). 

Here the Uj are given by the unique Cayley-Bacharach relation on A and l a is an extreme 
ray o/E| 6 . 

Note that the Cayley-Bacharach coefficients Uj can be computed by solving a system of 
linear equations (see exemplarily in Section [5]). Additionally, all Uj are rational, if every 
point in A is rational. Note furthermore that, for an arbitrary p, it is not clear whether 
an A with V(p) C A satisfying ( 13.1 j) does always exist. We discuss certain special cases 
in the two following sections. 

Proof. Let p G dPz$ with V{p) = {v%, . . . , V/-}, 1 < k < 8 such that the genericity condition 
( 12. 2p holds for V{p). We choose points V/~+i, ■ ■ ■ ,v$ such that ( 12. 2 p is still satisfied. We 
obtain Vg as the intersection of two relatively prime, cubic polynomials spanning up the 
vector space of all ternary cubics vanishing on V\, . . . , Vs (see Lemma [273]) . Notice that we 
obtain v 9 {vi, . . . , t> 8 } generically due to Lemma 12^1 and v 9 has to be real since Vi, . . . ,v 9 
is the intersection of two real polynomials (see e.g. [21]). Let u±, . . . ,Ug be the unique 
Cayley-Bacharach coefficients for v\, . . . , v s in the sense of (12. ip . Since t>i, . . . , v 9 G M 3 we 
have ux,...,Ug el (see [3 Lemma 4.1]). 

By Theorem 12 . 71 every vector a := (a±, . . . , a 9 ) el 9 satisfying (12. 5p with ai, . . . , a§ > 0, 



a 9 < yields an extreme ray l a : — > R, / S*=i a jf( v j) °f the dual SOS cone £3 



J 3,6' 



l a is the dual of a separating hyperplane for p if l a (p) < 0, i.e., since V(p) = {v\, . . . , v k }, 
if a k +ip(v + • • • + a 9 p{y 9 ) < 0. By ( 12. 5 p this is equivalent to 



a k+l p(v k+l ) + — 2 9 - -rp{vg) < 

1 [_ _| 8 



U l 1 1 M 8 \ ^ „.2„ 



(a fc+ ip(?; fe+ i) H h a 8 p(v 8 )) ■ 1 1 < u 9 p(v 9 ). 

V a i a sJ 

2 

Let \ ai ,...,a k ■= Ylj=i ir~ ' ( a k+iP(vh+i) H h 08P(fs)) > 0. Thus, Z a (p) < 0, if 



8 



i=fc+i \ ie{fc+i,...,8}\{i} 1 ) 

We choose a fc+1 := 1, . . . , a 8 := 1 and obtain 

\a u ...,a k + (U 2 k+1 + ■ ■ ■ + U 2 8 )(p(v k+1 ) + ■ ■ ■ + p(v 8 )) < ulp(v 9 ). 

Since lim aii ... iafc _ )>00 X ai ,...,a k \ the relaxation (13. ip yields an extreme ray l a on A sepa- 
rating p from the SOS cone S 3i6 . 

If on the other hand any SOS polynomial g G £3,6 with V(g) = V(p) would satisfy 
( 13. ip . then it follows from the upper construction that l a (g) < for a\, . . . , a k sufficiently 
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large, a k +i, . . . , a$ = 1 and a$ given by (12. 5p . This is a contradiction with Blekherman's 
Theorems 12.61 and 12.71 Thus, (13. ip indeed is a certificate for p to be not SOS. □ 

Note that it is also easy to show that condition (13.11) is never satisfied for SOS polyno- 
mials g G S 3 6 with V(g) = V(p) in order to prove extremality of l a by direct calculation 
and without usage of Blekherman's Theorems. This follows already from g = YTj=i ^1 
and the Cayley-Bacharach relations. 

In order to prove an analogon of the upper theorem for P^ \ £4,4, we need to show 
that Lemma [2.31 also holds for a seven point set A := . . . , t> 7 } C R 4 . Generically, the 
vector space of all quadrics vanishing on A has dimension three (see 

Lemma 3.2. Suppose A := {vi, . . . ,Vj} is a set of seven distinct points in R 4 , no four 
on a plane such that qi,q2,q3 is a basis for the vector space of all homogeneous quadrics 
with projective variety affinely represented by A. Then q±,q2,q3 are relatively prime. 

Proof. Suppose q±,q2,q3 have a common factor g. Then qj = g ■ q'- for j G {1,2,3} and 
g,q'j have to be linear in R[xi, . . . , x 4 ]. Due to the genericity condition at most three 
zeros (w.l.o.g. vi,V2,vs) are located on V(g) since otherwise there would exist at least five 
points are contained in a plane. Hence, V(q[), V(q' 2 ) and V(g 3 ) share four points which is 
a contradiction since for each j all points in V(q'j) are contained in a line. □ 

Theorem 3.3. Let p G dP^ such that p is not SOS . Let A := {v\, . . . , f§} C R 4 such 
that the genericity condition (12. 4p holds and V(p) = {t>i, . . . , v k } with 1 < k < 7 (i.e.: 
the projectivization of A is the intersection of three coprime polynomials qi, q 2 , <?3 G H^). 
Then there exists a certificate l a : if 4j4 — > M, / h> J2j=i a jf( v j)> a := ( a i5 • • • 5 a s) £ ^ 8 
with respect to A for p to be not SOS , if the following inequality holds 

(3.2) (ul +1 + --- + u 2 7 )(p(v k+1 ) + ---+ P (v r )) < u 2 8 p{v 8 ). 

Here the Uj are given by the unique Cayley-Bacharach relation on A and l a is an extreme 
ray o/£| >4 . 

The proof works the same way as for Theorem 13.11 with the obvious modifications. 
In fact, the proof of Theorem 13.11 already shows one possible way how to choose a = 
(di, . . . , ag) G M 9 to obtain a separating extreme ray l a . 

Corollary 3.4. For p G <9P 3)6 \ £3,6 and A = {vi, . . . ,t> 9 } D V(p) with 03.11) satisfied, 
one valid certificate is given by ai = ■ ■ ■ = a k = N e M. (for N sufficiently large), a k +i = 
• • ■ = as = 1 and a 9 given by (I2.5p . For p G (9P 4)4 \ £ 4>4 and A = {v\, . . . , v$} D V(p) with 
(13.21) satisfied, one valid certificate is given by a\ = ■ ■ ■ = a k = N G R (for N sufficiently 
large), a k +i = ■ ■ ■ = a 7 = 1 and as given by analogon of (12.51) for £44 (see \§\). 

In particular, l a is a rational certificate, i. e. every aj is rational, if every point Vj G A 
is rational. □ 

If one is interested in computing rational certificates, then there is the following problem 
from an application point of view. Suppose, we have a rational variety V(p) = {vi, . . . , v k } 
and we choose v k +i, ■ ■ ■ , v$ G Q 3 (resp. v k +i, ■ ■ ■ , v 7 G Q 4 ) such that the genericity condi- 
tion (12.21) (resp. (12. 4p ) holds (which is always possible). Then it is not clear a priori that 
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the ninth intersection point Vg G M 3 (resp. eighth intersection point Vs G M 4 ), given by 
Bezout, is rational, too. 

According to a recent preprint by Q. Ren [19], for p G P3,6 \ £3,6, the ninth intersection 
point can always be computed exactly. In particular, it can be followed that Vg will always 
be rational whenever v\, . . . ,v& are rational and hence, whenever V(p) is rational. 

Corollary 3.5. Let p G dP 3fi \ £ 3j6 with V(p) = {v 1 , . . . , v k } C Q 3 and {v k+1 , ...,v 8 } C 
Q 3 such that (12.21) holds. Then there is a rational certificate l a on A = {vi, . . . ,Vg} with 
vg given by Bezout, whenever (13. ip holds. □ 

Note that in our Theorems 13. II and 13.31 we only consider real points v 1, . . . , Vg whereas in 
Blekherman's Theorem 12.61 (at most) one pair of complex conjugated points are allowed. 
Anyhow, Blekherman conjectures in [5] that all extreme rays of the dual cones £3 6 and 
£4 4 can be described by nine (resp. eight) real points. Based on our results, we formulate 
a slightly different conjecture here. 

Conjecture 3.6. For p G P 3) e \ S 3i6 there exist Vi,...,Vg G M. 3 yielding a separating 
extreme ray forp in the sense of Theorem \2.6l Analogously, forp G P44 \ £4,4. 

Clearly, Blekhermans conjecture implies ours, since if every extreme ray is real rep- 
resentable, every nonnegative polynomial that is not SOS can be separated by a real 
intersection. But it is unclear whether the two variants are equivalent or not and, if not, 
whether our conjecture is a strong relaxation of Blekherman's. 

4. The seven point case 

Let p G Ps t Q \ £3,6 and assume we are interested in finding a separating extreme ray l a in 
£3 6 for p. That means we need to find a generic 9-point set A being the intersection of the 
varieties of two coprime polynomials qi,q% G H 33 such that the conditions in Theorem 
12.61 hold. If p is located on the boundary of P 3j6 , i.e. V{p) = {v 1 , . . . , v k } ^ 0, then 
Theorem 13.11 and Corollary 13.41 yield a certificate l a whenever one can fill up V(p) with 
points v fc+i, . . . , Vg such that condition (13. ip is satisfied. 

Anyhow, it is not obvious a priori if resp. how points v k +i, ■ ■ ■ ,vg can be chosen such 
that the sufficient condition (13.11) of Theorem 13.11 holds (note that vg is always given by 
Bezout in this approach). It turns out that for k = 7, i.e. the easiest non-trivial case, 
the problem of choosing an appropriate v$ is easy since almost every v$ yields a v 9 such 
that (13. ip is satisfied. Note that for p G dP^^ \ S 3 6 the variety V(p) always satisfies the 
genericity condition (12. 2p (see [2Tj). 

One reason why this case might be of special interest is that k = 7 with v\, . . . ,v>? 
satisfying the genericity condition (12. 2p is the smallest number of zeros of a nonnegative 
polynomial such that the dimensional difference between exposed faces of P 3 6 and S 3 6 
given by the vanishing of polynomials on these zeros is strictly positive (see jl]). 

Theorem 14. II holds analogously in the case of p G P44 \ £4,4 and k = 6 with the obvious 
modifications. We omit to formulate the result for this case separately. 

Theorem 4.1. Let p G dP^^ \ £3,6- Let A := {vi, . . . ,Vg} C M 3 such that the genericity 
condition (12.21) holds and V(p) = {t>i, . . . ,^7}. Then there exists a certificate l a : H 36 — > 
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K, / i— > Y^j=i a jf( v j)> a j ^ with respect to A for p to be not SOS if 
(4.1) u 2 8 p(v 8 ) ^ ulp(v 9 ). 

Furthermore, l a is an extreme ray o/E^g. 

Proof. We choose v% such that the genericity condition ( 12. 2ft still holds for {vx, . . . ,Vg} 
and obtain a real v 9 {v%, . . . ,v$} such that ( 12. 2 p and Cayley-Bacharach holds for A 
generically (see proof of Theorem 13. ip . Since all conditions of Theorem 13.11 are satisfied, 
there is a certificate I if u1p{v&) < Ugp(vg). But since {v±, . . . , i>7,i>g} also yields A with 
Bezout's theorem, this condition holds w.l.o.g. as long as u 8 p{y%) ^ Ugp{vg). □ 

Corollary 4.2. Conjecture \3.6\ holds for p G dP^e \ £3,6 with #V(p) > 7 and p G 
aP 4 ,4\£ 4) 4 wiih#V(p) > 6. 

Let p G P3 t 6 such that p is not S'OS' with V(p) = . . . ,Vk} C M 3 . Assume that we 
want to test for a large amount of point sets {v k+i, ■ ■ ■ , v$} in R 3 if (13. ip can be satisfied 
for each point set (under the condition that (12. 2p holds for v\, . . . ,t>s). We have to do 
two steps. Firstly, we have to compute vg by calculating two coprime polynomials qi,q2 
vanishing on {v\, . . . , Vs} via solving a system of linear equations on the coefficients of a 
ternary cubic and then calculate their Grobner basis. If we are in the case (n, 2d) = (3, 6), 
then this can also be done by using the formula given in [15] . Secondly, we have to 
compute the Cayley-Bacharach relation on A := {vi, . . . ,Vg} via solving a system of 
linear equations again. 

In the seven point case, i.e. if we want to check ( 14. ip for a large amount of point sets, 
we can avoid calculating the Cayley-Bacharach relation for each {v$,vg} if we are not 
interested in computing the separating hyperplane given by l a explicitly. The trick is to 
use Reznick's Lemma 1231 Let V(V(p)) := {h G -£^3 : V{p) C V(h)} be the vector space 
of all real ternary cubics vanishing on V(p). By Lemma 12.51 there is a basis {hi, /i2, ^3} 
of V(V{p)) such that V{h x ) H V{h 2 ) H V{h 3 ) = V(p) and hi, h 2 , h do not depend on the 
choice of v 8 and t> 9 . 

Corollary 4.3. Let p and A be as in Theorem \4-l\ Let h G V(V(p)) with V(h) = V(p). 
Then there exists a certificate l a : H 3 6 — y R, / !->■ Ylj=i a jfi v j)> a j e ^ with respect to A 
for p to be not SOS if 

h 2 (Vg)p(v 8 ) ^ h 2 (v 8 )p(Vg). 

Furthermore, l a is an extreme ray o/Sgg. 

Proof. By Lemma 12.51 there is a h G V(V(p)) with V(h) = V(p). Since A satisfies (12. 2p 
and h G 3 the Cayley-Bacharach relation ( 12. ip holds for h. Since V(h) = V(p) = 
{v 1, . . . , vj} this means u&h(v%) + ugh{vg) = and thus u\h 2 {ys) = u1h 2 {yg). Hence, 

ulp(v 8 ) ^ U 2 9 p(vg) <^> h 2 (vg)p(v 8 ) ^ h 2 (v 8 )p(vg). 

The assertion follows with Theorem 14.11 □ 

If we want to compute a specific l a the approach of Corollary 14.31 does not suffice since 
by (12. 5p we need to know the coefficients Uj of the Cayley-Bacharach relation (12. ip on A 



SEPARATING INEQUALITIES FOR NONNEGATIVE POLYNOMIALS THAT ARE NOT SOS 11 



to compute the scalar ag for l a . 

In the appendix we present an example of our method applied to a nonnegative poly- 
nomial with exactly seven zeros that is not SOS. 

5. AN APPLICATION: THE MOTZKIN POLYNOMIAL 

In this section we demonstrate applications of our method. It turns out that finding 
a separating extreme ray becomes more difficult for polynomials in p G dPs^ \ £3,6 with 
six or less zeros. The condition (13. ip provides more degrees of freedom and, in particular, 
the left hand side of this inequality has more than one summand. This fact yields that 
the set of point configurations A := {vi, . . . , vg} C M 3 with V(p) C A which do not satisfy 
(13. ip will not be a lower dimensional subset in general - in contrast to the seven point 
case (independently from which point corresponds to the negative entry ag in an extreme 
ray). The same difficulties arise for p G P^4 \ £4,4 with five or less zeros. 

Exemplarily, we investigate the Motzkin polynomial 

/ \ 42,24 o 2 2 2 , 6 

m{x } y } z) = xy+xy— Sxyz+z. 

The Motzkin polynomial has six zeros 

Vl := (1,0,0), v 2 := (0,1,0), v 3 := (1,1,1), 
v 4 := (-1,1,1), v 5 := (1,-1,1), v 6 := (1,1,-1). 

As a first instance we choose t> 7 := (0, 4, 1) and v 8 := (4, 0, 1). These eight points satisfy 
the genericity condition (12. 2p . This can be checked by looking at the 3x3 minors of the 
8x3 matrix given by the coordinates (x, y, z) of the points v\,...,v% and looking at the 
6x6 minors of the 8x6 matrix given by the coordinates (x 2 ,y 2 , z 2 ,xy,xz,yz) of the 
points v 1, . . . , Vs. Hence we can compute two coprime ternary cubics 

o 9 9 9 9 

qi = —16z + Ihx z + y z + 56y x — 56z x 

00 q 9 9 

q 2 = —Az — x y + Ax z + loy x — VSz x + z y 
vanishing on Vi, . . . , Vs by solving the system of linear equations 

h{vi) = 0,...,h(v 8 ) = 
on the coefficients of h, where 

h := bix 3 + b 2 y 3 + b 3 z 3 + b^x 2 y + b 5 x 2 z + b 6 y 2 z + b 7 y 2 x + b 8 z 2 x + b 9 z 2 y + b w xyz 
again. We compute the Grobner basis 

{7z - 26z 2 - 15z 3 + 26z A + 8z 5 , -2z + 8z 2 + 2z 3 - W5y - 8z A + W5z 2 y, 
8z 4 - A22z 3 + 1575y 2 - 1583z 2 + 422z, x - 1} 

of qi, q 2 and x — 1 with lexicographic ordering. We obtain t> 9 = (1, 1, —7/2) and compute 
the Cayley-Bacharach coefficients Uj by solving the system of linear equations 

U\h{vx) + ■ ■ • + u 9 h(vg) = 
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in «i, • • • , ug. The solution is (up to multiplicity) 

( ,a 40 . „ 24 1 1 H8098\ T 

u = _64,-64, ,-4, -4, —,1,1, 

V 9 ' ' ' 5 ' ' ' 5 J 

We have 

(it? + Wg) ■ (m(f 7 ) + m(v 8 )) = 4, 

w 2 m(t; 9 ) = 228. 

Hence, the condition (13 .ip of Theorem 13.11 is satisfied and we find a separating hyper- 
plane for m on A. According to Corollary 13.41 we choose a\, . . . , ciq := 100 and a-?, a§ := 1. 
By (12. 5 p we obtain 

— Ug -14121476824050 
ag ~ d + + ~ 2143157 ' 

We check the correctness of our result by 

la{m) = axm{v\) + ■ • • + agm(vg) 

N 14121476824050 / 1 7 

= m(0,4, 1) + m(4,0,l) m 1,1, — 

v ' ' ; v ' ' J 2143157 V 2 

1484936 

< 0. 



2143157 



Thus, by Blekherman's Theorem 12. 7\ l a is a (rational) extreme ray of Sg 6 separating the 
Motzkin polynomial m from the SOS cone. 

In contrast to the seven point case, not every generic point configuration yields a sep- 
arating certificate. For example, with the same approach it is easy to show that the 
instance v-j := (2/7, 2/3, 1) and v$ := (2/3, 2/7, 1) does not satisfy the condition (13.11) (see 
Appendix IBl. 

We show that for a symmetric choice of vj,vs, i.e. v-j = (q,s,l),v$ = (s,q,l) with 
g,sel, the set 

S := {{q,s)eR 2 : (v% + u 2 8 ) ■ (m(v 7 ) + m(v s )) < u 2 g m{v 9 )} 

yieding a 9-point configuration which satisfies (13 .ip is full dimensional with some nice 
geometric structure (see Figure (JTJ)). 

With the formula in [19] we obtain the ninth Cayley-Bacharach point 



7 3 s + 2q 2 s 2 — q 2 + s 3 q — 2sq - s 2 
u 9 = 177 \ 'I ^ 3fi + 2g ,2 s 2 - q 2 + s 3 q - 2sq - s 2 

„3 I „2„ _ n„ i „2„ _ o„ i „3 



+ O - 2g + S^Q - 2s + S d 



n{q,s) 
where 

n(q, s) = 12 • (g 3 s 2 - 2q 3 s + q 3 + s 3 q 2 - 2q 2 s 2 + q 2 s - 3q + 4sg - 2s 3 g + s 2 q + 2 - 3s + s 3 ) 
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vanishing at q = —2 — s,q = l,s = 1. Furthermore we obtain a non-generic point set 
for q = s,q = —s,q = —l,s = —1 and q = 2 — s. We compute the Cayley-Bacharach 
coefficients in dependence of q, s and obtain 



Note that the numerator of (it 2 . + u|) • (m(t> 7 ) + m(t> 8 )) is exactly the dehomogenized 
Motzkin polynomial in s,q for z — 1. We set 



Since g 7^ s by assumption we just need to investigate K(q, s). Thus, by (13.1 p we have 
(q, s) G S if and only if if(g, s) < and q £ {±s, ±1, ±2 — s}. Equivalently, S 1 = if and 
only if K(q, s) is PSD and q ^ {±s, ±1, ±2 — s}. Since s) is a bivariate polynomial 
of degree 4, if (g, s) is PSD if only if it is SOS. It can be checked easily that this is not 
the case. 

We provide a plot of the set S := {(q, s) : K(q, s) < 0}\{(q,s) : q = 2 — s or q = — 2 — s} 
in Figure [1] (note that the other non-generic cases are not part of S although they are 
relevant for the computation). Obviously, this set is symmetric in q = s and q = —s, 
semi-algebraic and for every q there is an s such that (q, s) G S. 



(uj + it|) ■ (m(vr) + m(v 8 )) 



UQm(vg) 



64(1 + s A g 2 + g 4 s 2 - 3gV) 
(q - a y 

16(2<j 4 s 2 + q 4 + 4s 3 q s - 4q 3 s + 2s 4 q 2 - 6q 2 s 2 - 2q 2 - 4s 3 q + 4sq + s 4 - 2s 2 + 4) 



K(q,s)/L(q,s) 



{u 2 7 + u\) ■ (m(v 7 ) + m(v&)) - ulm(vg) 
16(2g 2 s 2 - q 2 + 2sq - s 2 + 2)/(q - s) 2 . 




FIGURE 1. The set S := {(q, s) : K(q, s) < 0} \ {(q, s) : q = 2 - s or q 
—2 — s} is given by the red area without the blue lines. 



Due to the rich geometric structure of S, it would be of interest to investigate the geo- 
metric structure of the set of appropriate point configurations satisfying ( 13. ip for general 
nonnegative polynomials with k zeros. 
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In contrast, we briefly demonstrate the numerical method for finding a 9-point certificate 
given in Section [Z2l and the corresponding problems. Let r = (x 2 + y 2 + z 2 ) 3 G int(£ 3) 6). 
Then, consider the following semidefinite optimization problem 

min A such that m + \(x 2 + y 2 + z 2 ) 3 G £ 3 6 . 

Agffi 

The optimal lambda is given numerically by A ~ 0.004596411406567 and the correspond- 
ing sum of squares decomposition of m + Ar « s\ + s 2 . + s 2 is given by 

m + X(x 2 + y 2 + z 2 ) 3 » (-0.858558xz 2 + 0.941354xy 2 + 0.067796x 3 ) 2 + 
(-0.858558?/z 2 + 0.067796y 3 + 0.941354x 2 y) 2 + (-1.002295z 3 + 0.360838y 2 z + 0.360838x 2 z) 2 . 

Now, one has to choose an appropriate linear combination of two of the polynomials s$ 
to obtain the nine intersection points. But it seems unclear how to do this. E.g., since 
the coefficients of the Si are given numerically, computing the Grobner basis of two of the 
Si and, say, x — 1 cannot be expected to work properly. 

Finally, we remark that our method allows to generate strictly positive polynomials that 
are not sums of squares including a certificate without optimization, if (13. ip is satisfied for 
a polynomial p G <9P 3]6 \ E 3)6 (resp. for (n, 2d) = (4,4)). Let l a be a separating extreme 
ray for p and define n := p + A ■ r with A G M>o and r G int(E 3)6 ) (resp. r G int(£ 4)4 )). 
Then n is strictly positive and evaluating n on l a yields 

la(n) = l a (p) + A • l a (r) 

with l a {p) < and l a (r) > 0. Hence, we can immediately solve for A such that l a {n) < 0. 
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Appendix A. An example for the seven point case 

As an example of our method we investigate the polynomial 

p := x 2 y 2 {x — y) 2 + y 2 z 2 (y — z) 2 + z 2 x 2 {z — x) 2 + xyz{x — y)(y — z)(z — x), 

given by Reznick in [21] . It has the seven zeros 

vi := (1,0,0), v 2 := (0,1,0), v 3 := (0,0,1), v 4 := (1,1,0) 
v 5 := (1,0,1), v 6 := (0,1,1), v 7 := (1,1,1). 

We choose v$ : = (—2, 5,-1) which satisfies (12. 2(1 . 

h := bix 3 + b 2 y 3 + b 3 z 3 + b±x 2 y + b 5 x 2 z + b 6 y 2 z + b 7 y 2 x + b 8 z 2 x + b 9 z 2 y + b w xyz. 

We compute two ternary cubics 

qi = —x 2 y — 35x 2 z + y 2 x + 35z 2 x, q 2 = I5x 2 z — y 2 z — 15z 2 x + z 2 y 

vanishing on v\ , . . . , v% by solving the system of linear equations 

h{ Vl ) = 0,...,h(v 8 ) = 

on the coefficients of h. qi and q 2 are coprime since the genericity condition ( 12 .2(1 holds 
for v\, . . . ,v$ by Lemma 1231 We obtain a ninth intersection point Vg = (3, 10, 1) of qi, q 2 
via computing the Grobner basis 

{6z 4 - llz 3 + 6z 2 - z, z 2 y + 35z 3 - zy - 50z 2 + 15z, -y - 35z + y 2 + 35z 2 , z - 1} 
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of qi, q2 and z — 1 with respect to lexicographic ordering. To compute the Uj of the 
Cayley-Bacharach relation (12. ip we solve the system of linear equations given by 

U\h{vi) + ■ ■ ■ + Ugh(vg) = 0. 

We obtain u = (84, -1260, -36, -90, 63, 35, -60, 3, 1) T and therefore 

u 2 8 p{v 8 ) = 48456, u 2 9 p{v 9 ) = 56016, 

i.e. (14. ip holds and we find a separating hyperplane for p on A. According to Corollary 
13.41 we choose oi, . . . , 07 := 10 9 and as := 1. By (12. 5p we obtain 

-ul -10 9 

119 ~ + ~ 84 2 + (-1260) 2 + (-36) 2 + (-90) 2 + 63 2 + 35 2 + (-60) 2 + 9 ■ 10 9 

ai ' 'as 

500000000 



4500806423 
We check correctness of our result by 

, / x /x / \ / \ 500000000 

la[P) = aip(t>i) + • • - + a 9 p{v 9 ) = p{-2, 5, -1) - 4500806423 • P( 3 > 10 > !) 

471957277321 

< 0. 



5626008028750 

Hence, by Blekherman's Theorem 12.71 Z a is an extreme ray of S3 6 separating p from the 
SOS cone. Since it is constructed in the way of Corollary 13.51 this certificate is indeed 
rational. 



Appendix B. A non- appropriate point configuration for the Motzkin 

POLYNOMIAL 

We choose the six zeros of the Motzkin polynomial and the instance t> 7 := (2/7, 2/3, 1) 
and v 8 := (2/3,2/7,1). These eight points satisfy the genericity condition (12.2]) . We 
proceed the same way as before and obtain two coprime ternary cubics 

qi = 272z 3 + 520x 2 y - A61x 2 z + 189y 2 z - 520z 2 y 

q 2 = -3A0z 3 - 4Qlx 2 y + 340x 2 ^ - 189y 2 x + 189z 2 x + 4Qlz 2 y 

vanishing on v\, . . . , i>g. The Grobner basis for qi, q2 and x—1 with respect to lexicographic 
ordering is 

{1365z - 2014z 2 _ 425z 3 + 1878z 4 - 940z 5 + 136z 6 , -329154z + 362620z 2 + 266458z 3 
-257985y - 362620z 4 + 62696z 5 + 257985z 2 y, -239454894z + 118408655z 2 + 210552038z 3 
-167167820/ + 28902856z 5 + 48759165y 2 , x - 1} 

and the ninth point is given by v 9 = (1,1, 65/34). The Cayley-Bacharach coefficients are 

/-264 -264 891 -99 -99 43659 43659 -4913 \ T 
U V 7 ' 7 'IT'TF'IF' ' 160 ' 160 ' 62 J 
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But now 

/ 2 2n / / x /xx n 1906108281 n 177025 n „ nn „, Ml 

(ul + ul) ■ (m(v 7 ) + m w 8 = 2 2 = 271097.1914 

v 7 v v t; v »jj 25600 19448I 

i*i-m(u 9 ) = 250270.0664. 

If we exchange the points v 7 and t> 9 or the points v$ and v 9 we obtain in both cases 

(i4 + u 2 8 ) ■ (m(v 7 ) + m(v 8 )) = 3291366.873 and u\ ■ m{v 9 ) = 67774.29785. 
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